Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 19 February 2008 



(MN style file v2.2) 



Daily Observations of Interstellar Scintillation in 
PSR B0329+54 



N. Wang,i* Z. Yan,^'^ R. N. Manchester,^ and H. X. Wang,^ 

^ Urumqi Observatory, NAOC, 40-5 South Beijing Road, Urumqi, Xinjiang, China, 830011 
'^Australia Telescope National Facility, CSIRO, PO Box 76, Epping, NSW 1710, Australia 
^Graduate University of Chinese Academy of Sciences, 19A Yuquan road, Beijing, China, 100049 



19 February 2008 



OO 

o 
o 

(N 

CD 

S: 

9^ 
6 

> 

OO 

m 
\6 
o 
i> 
o 

> 
X 

1 INTRODUCTION cussing effects observed as RISS. A goo d rev iew of the theory 

■ " " ' behind this interpretation is given by ,Ricketti (il99Qi ) . 

In the early days of pulsar observations, slow fluctuations • n . i .i . .i . r i . 

r— — -1 I -] It IS generally accepted that the spectrum ol electron 

of pulsar amplitudes were reported by ILole et al.l (11970 ) j a ^ i- ■ • ^ ^ n j- i 

f rr; 71 J 71. ^ — '■, — density nuctuations m the interstellar medium has a power- 

and lHuguenin et al.l l| 19731 ). But it was not recognized that form- 

they were an interstellar propagation effect until a close 

correlation between the time-scale of the intensity fluctua- Psniq) = C^q'^ , (1) 

tions and the pulsar dispersi on measure was poin ted out by 

[Sicb or (1982). Subsequentlv. [Rickett et al.| (jl984| ) explained where q = 2'k/L is the wave- vector associated with a spatial 
these slow fluctuations as refractive interstellar scintillations size L. The quantity Cn is a measure of turbulence along a 
(RISS). This phenomenon is produced by the same spectrum particular hue of sight. The power-law index f3 which indi- 
of electron density fluctuations in the interstellar medium cates the steepness of the inhomogeneity spectrum is in the 
that is responsible for the weU known diffractive scintilla- range of 3 < /3 < 5. But the details of the spectrum, e.g. 
tion (DISS). DISS is caused by the small spatial scale den- its slope and the range of scale sizes over which it is valid, 
sity fluctuations (10® — 10* m) and it appears as intensity stiU remain open questions. Also, different assumptions are 
variations in both the time and frequency domains with made about the distribution of scattering material along the 
characteristic scales of minutes to hours and kHz to MHz path to the pulsar. Many observations are well accounted for 
respectively. On the other hand, large-scale inhomogeneities by models in which the scattering is concentrated in a single 
(10^° — 10^^ m) in the interstellar medium give rise to fo- centrally located thin screen, but others imply that there 

are many scattering sites or a statistically uniform screen 
extending over most or all of the path. 

A range of observations suggest that a Kolmogorov fluc- 
* Emailina. wang@uao.ac.cn tuation spectrum with (3 = 11/3 is widely applicable in the 



ABSTRACT 

Quasi-continuous observations of PSR B0329-I-54 over 20 days using the Nanshan 25-m 
telescope at 1540 MHz have been used to study the effects of refractive scintillation on 
the pulsar flux density and diffractive scintillation properties. Dynamic spectra were 
obtained from datasets of 90 min duration and diffractive parameters derived from 
a two-dimensional auto-correlation analysis. Secondary spectra were also computed 
but these showed no significant evidence for arc structure. Cross correlations between 
variations in the derived parameters were much lower than predicted by thin screen 
models and in one case was of opposite sign to the prediction. Observed modulation 
indices were larger than predicted by thin screen models with a Kolmogorov fluctuation 
spectrum. Structure functions were computed for the flux density, diffractive timescalc 
and decorrelation bandwidth. These indicated a refractive timescale of 8 ± 2 h, much 
shorter than predicted by the thin screen model. The measured structure-function 
slope of 0.4 ±0.2 is also inconsistent with scattering by a single thin screen for which a 
slope of 2.0 is expected. All observations are consistent with scattering by an extended 
medium having a Kolmogorov fluctuation spectrum which is concentrated towards 
the pulsar. This interpretation is also consistent with recent observations of multiple 
diffuse scintillation arcs for this pulsar. 

Key v^fords: pulsars:individual:PSR B0329+54 - ISMrstructure 
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interstellar medium (| Armstrong et ai1ll995h . However, ob- 
servational studies showing, for example, persistent drift- 
ing bands in dynamic spectra and increased modulation of 
diffractive properties suggest that excess power is often seen 
at low spatial frequencies, implying a spectral index /3 4. 
Exam ples of extensive studies are those by I Gupta et al.l 
1 19941 ) ■ who studied eight pulsars at 408 MHz using the 
Lovell telescope and interpreted changing diffractive pat- 
terns in dynamic spectra in terms of refractive effe cts, and 
those bv lBhat et all (Il999al lbh and lBhat etHI (|l999l ). These 
latter authors used the Ooty radio telescope at 327 MHz to 
study 18 pulsars over a 2.5-year dataspan and showed that, 
while many observations were consistent with a Kolmogorov 
spectrum, other s show ed evidence of a steeper spectrum. 
IStinebring et al.1 (|2000l ) reported on 5 years of daily moni- 
toring of the flux density of 21 pulsars at 610 MHz, showing 
that most of the results, covering a wide range of scattering 
strengths, were consistent with a Kolmogorov fluctuation 
spectrum. In some pulsars, though modulation indices were 
higher than expected, possibly due to the presence of an 
"inner scale" , a cutoff in the fluctuation spectrum at scales 
greater than the diffractive scale Sd ~ A/6d, where A is the 
wavelength and 6d is the angular size of the scattering disk. 

Other observations suggesting steeper power-law 
spectra (/3 ^ 4) include fringing in d ynamic spec- 
tra ir nplying multiple ir naging (e.g., [ Cordes et al.l 



1 19861 : ICordes fc WolszczanI Il986l: iRickett et al 



and extreme scatter i ng e vents (jLestrade et al 



.19971 ) 

.. Il998h . 

iBlandford &: NaravanI l|l985l ) suggested a steeper spectrum 
to deal with the theoretical difficulties in supporting a 
turbulent cascade in the Kolmogoro v spectrum. Ba s ed on 
power-law models of the spectrum, iRomani et al.l (|l986l ) 
analysed the effects of RISS on diffractive parameters 
for different slopes /3 = 11/3, 4 and 4.3, predicting anti- 
correlations for (_F,A!-'d) and [F, Aid), where F is flux 
density and Ai/d, Aid are the diffractive scintillation band- 
width and timescale, respectively and a postive correlation 
for (Ai'djAfd), with the magnitudes of the correlation 
coefficients being larger for steeper spectra. 

More recent observations have tended t o concentrate 
on de t ailed studies of individual pulsars fe.g.. | Shishov et al] 
I2OO3I : iRamachandran et all bood : ISmirnova et al.1 12OO6I ) or 
on the fascinating "scintillation arcs" which are seen 
in secondary spectra, that is, two-dimensional Fourier 
transforms of dynamic spectra, in high-sensit i vity, high- 
resolution observation s (jStinebring et all 1200 ll : iHill et al] 
I2OO3I : [ Stinebrind 1200^ ). These arcs, which are closely re- 
lated to multiple imaging and the frequently observed 
"criss-cross" sloping bands, result from interference between 
rays in a central core and rays from an extended scat- 
tering disk. They form a powerful probe of structure in 
the interstellar medium (jCordes et al.ll2006l ). For example, 
IPutnev fc Stinebrind (|2006l ) present recent observations of 
six pulsars showing multiple arcs of different curvature im- 
plying a distribution of scattering centres along the line of 
sight. 

PSR B0329-I-54 is one of the strongest pulsars known 
and its scintillation properties have been studied by many 
autho rs. Observations at 610 MHz by IStinebring et al] 
(|l996t ) showed that the correlations between variations of 
flux, decorrelation bandwidth and scintillation time-scale 
were consistent with the theoretical predictions with /3 < 4. 



However. [Bhat et al.1 l|l999bl ) found that for PSR B0329+54, 
correlations between these parameters at 327 MHz were not 
in accord with the predi ctions. Long-term flux density mon- 
itoring at 610 MHz bv iStinebring et al.l (|2000l ) showed that 
the modulation index was somewhat higher than expected 
for a Kolmogorov spectrum, i.e. 13 > 11/3. IShishov et al.l 
l|2003l ) took data obtained over a wide range of frequencies 
flnding a value of f3 — 3.5, less than but consistent with a 
Kolmogorov turbulence spectrum, and interpreted this re- 
sult in terms of weak plasma turbulence in the interstellar 
mediu m. Observations a t 1540 MHz using the Nanshan tele- 
scope (|Wang et al.ll2005l ) showed inconsistencies in values of 
the power-law index obtained from two different approaches. 
The ratio of refractive scattering angle 6r and diffractive an- 
gle 6d gave /3 « 3.6, while the frequency dependence of scin- 
tillation parameters gave /3 > 4. These conflicting results 
suggested undertaking observations with a better sampling 
of the refractive variations. These have a predicted timescale 



Atr ~ {2iy/Ai'd)Atd, 



(2) 



(|Gupta et al.lll993l : [Stinebring et al.ll2000l ). about four days 
for this pulsar. 

In this paper, we present the results of quasi-continuous 
1540 MHz observations of PSR B0329-I-54 over a 20-day in- 
terval in 2004, March. Diffractive parameters were sampled 
at 90-min intervals, giving good resolution of the expected 
RISS variations. The layout for the rest of this paper is as 
follows: in Section 2 we introduce the observations; data 
analysis and results are shown in Section 3, and in Section 
4, we discuss the results and their interpretation. Conclu- 
sions are presented in Section 5. 



2 OBSERVATIONS 

The observations were made using the Nanshan 25 m radio 
telescope with a central frequency of 1540 MHz. Our system 
is sensitive to orthogonal linear polarizations and the system 
temperature for the two channels is 20 K and 22 K respec- 
tively. Each polarization has 128 channels, each of band- 
width 2.5 MHz, giving a total bandwidth of 320 MHz. After 
detection and high-pass filtering, the data are 1-bit sampled 
and folded at the apparent pulsar period. Signals from the 
two polarizations are summed to form total intensity pulse 
profiles. 

We made our observations from March 12 to March 
31 in 2004, and obtained more than 150 scintillation 
dynamic spectra for PSR B0329-f54. The observations 
were continuous (PSR B0329+54 is circumpolar at Nan- 
shan) apart from an interruption of about 3 days from 
March 20 due to telescope problems. In our analysis, 
the time block for each dynamic spectrum is typically 
90 minutes with a sub-integration time of 60 s. We cal- 
ibrated the pulsar flux density scale using observations 
of ten strong, relatively distant pulsars which all have 
well-measured flux densities in the ATNF pulsar catalogue 
(www.atnf.csiro.au/research/pulsar/psrcat). 
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3 DATA ANALYSIS AND RESULTS 

3.1 Scintillation Dynamic Spectra 

Dynamic spectra were obtained for each observation by plot- 
ting the mean flux density in each channel against time. Ex- 
amples of these plots are shown in Fig. [T] Channels badly 
affected by radio-frequency interference have been interpo- 
lated over from adjacent frequency channels. 

Fig.[l]shows that the dynamic spectra are well resolved 
and that the scintles or peaks of flux density change their 
shapes and sizes in the frequency and time domains. At 
some epochs, the pulsar seems weaker and the scintles are 
smaller (MJD 53088.4), and at other epochs, the pulsar sig- 
nal is stronger and the scintles larger (e.g. MJDs 53081.6, 
53083.2, 53093.3). The frequency dependence of DISS can be 
seen within our receiver bandwidth, i.e., scintles are smaller 
at the lower end of the observed band (MJD 53091.8). At 
some epochs the spectra show a significant drifting pattern 
(e.g. MJD 53082.4), which implies a modulation by RISS. 
H owever, we d i dn't o bserve the fine fringes seen in the data 
of I Wang et ai] (|2005l ). The scintles are generally smaller in 
this observing session, resulting in a smaller value in Ai/^, 
which we will discuss more in Section r3.4l Deep modulations 
are common in all the dynamic spectra we have recorded. 

The limited number of scintles in the dynamic spectra 
introduces a statistical estimati on error in the scintillation 
parameters l|Cordes et al.|[l985h . An approximate estimate 
of the number of scintles N in the dynamic spectra is given 
by: 

_ Tobs X BWobs ,„s 

^ Aud X Aid ' 

where Tobs is the total observing time and -BW^oba the total 
observing bandwidth, Avd decorrelation bandwidth and Atd 
is the DISS timescale. The fractional estimation error is then 
given by: 

= " A.dxAU ) ^ 

where we assume a scintle filling factor of 0.5 (|Bhat et al.l 
This error contribution is taken into account in the 
estimates of scintillation parameters in Section [3.41 



3.2 Secondary Spectrum 

The two-dimensional Fourier spectrum of the primary dy- 
namic spectrum is often referred to as the seco ndary spec- 
trum . To show the beautiful arc structures (e.g.. Istinebrind 
I2OO6I ). high frequency resolution and high sensitivity are re- 
quired. Unfortunately, the Nanshan system does not pro- 
vide these. H owever, earlier obse rvations at 1540 MHz for 
this pulsar by IWang et all (|2005h showed sloping fringes in 
the dynamic spectra resulting in offset features in the sec- 
ondary spectra. Surprisingly, despite using the same system 
at Nanshan, our three-week consecutive observations show 
no significant arc structure or offset features in the secondary 
spectra. 



3.3 Two-Dimensional Auto-Correlation Function 

To quantify the diffractive parameters, we use the two- 
dimensional auto-correlation function (ACF), which was 



computed for frequency lags up to half of the observing 
bandwidth and for time lags up to half of the observing 
time. The ACF of the dynamic spectra A{Ai', At) is defined 
as: 

A{Au, ^t) = J2J2 ^-^('^' + Au,t + At), (5) 

u t 

where AF{i', t) = Fiy, t) — F, and F is the mean pulsar fiux 
density over each observation, then the normalized ACF is: 



p{Au, At) = A{Au, At)/A{0, 0). 



(6) 



We plot the normalized ACF and one dimensional cut 
at zero time and frequency lag in Fig. [51 These figures cor- 
respond to the dynamic spectra shown in Fig. [T] Following 
convention, the DISS time-scale Atd is defined as the time 
lag at zero frequency lag where the ACF is 1/e of the max- 
imum, and the decorrelation bandwidth Aud is defined as 
the half-width at half-maximum of the ACF a long the fre- 
quency lag axis at zero time lag (jCordeslll986l ). For weaker 
pulsars, random system noise results in a strong spike at the 
zero time and frequency lag of ACF. For PSR B0329-I-54 
the spike is weak but still observable at some epochs. To 
remove the spike at the origin we used the four neighbour- 
ing frequency-lag points to fit for a parabola across the zero 
frequency lag. The central point is interpolated from this fit. 



3.4 ACF and Scintillation Parameters 

Following ICupta et all (jl994l ) and iBhat et all l|l999al ). we 
use a two-dimensional elliptical Gaussian function to fit the 
ACF with form: 

p{Aiy, At) = Co exp[-(Ci Aiy^ + CiAvAt + C^At^)], (7) 

in which Co is unity since the ACF is normalized to unity. 
By using a minimization procedure, we obtained param- 
eters Ci, C2 and C3. The fitting procedure was a search for 
the three non-linear parameters for the least-squared error 
between model and data over the central region of ACF. The 
scintillation parameters Avd and Atd are calculated as: 



Ai^d = \/ln2/Ci, 
Atd 



(8) 
(9) 



Following iBhat et al.l (|l999al ). we use dt/dv to describe 
the orientation of the elliptical Gaussian. This is propor- 
tional to the refractive scattering angle 6^, and is given by: 



dt_ _ / C2 \ 
du ^ V2C3/ 



(10) 



From the analysis, we also obtained the uncertain- 
ties for Ci, C2 and C3. Using these uncertainties added in 
quadrature with the statistical error (Equation|3}, we calcu- 
late the errors of Avd, Atd and dt/dv. Since dt/du has zero 
mean, we applied the statistical error to (IC2I + \^CiCs) /2 
rather than to C2 itself. For flux density, measurement errors 
are estimated from the baseline noise in the mean pulse pro- 
files. Time variations of F, Aud, Atd and dt/dv along with 
their uncertainties are shown in Fig. [3] These parameters 
show significant variations with time which can be explained 
as an RISS effect as discussed in Section |4l 
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Figure 1. Samples of dynamic spectra for PSR B0329+54 taken from Marcli 12 to March 31 in 2004. The horizontal and vertical axis 
corresponds to frequency and time respectively. Darker regions of the grey-scale saturate at black for intensities above 85 percent of the 
maximum, white corresponds to below 1 percent of the maximum intensity and the grey-scale is linear between these limits. The Modified 
Julian Date (MJD) of the observation is shown in the top right corner of each panel. 



3.5 Derived Scattering Parameters 

Given a power-law spectrum of electron density fluctuations 
(Equation [1]), the average scattering level along the line-of- 
sight, Cn, for a Kolmogorov spectrum is given by 



eters is: 



Do 



-11/6 



-20/3 
(11) 



l|Cordes et al.lll985l ). 

Another quantitative measurement of the strength of 
scattering is the parameter u which is defined as the ratio 
of the Fresnel scale sf to the coherence scale Sd = («:&d)~ 
l|Rickettlll990l ). where 9^ is the diffractive scattering angle. 
For strong scattering, u ^ 1. In terms of the diffractive 
scintillation parameter Ai/^: 



2v 



(12) 



The diffractive timescale depends on the velocity of the 
ray path across the scattering screen. Assuming a thin screen 
and that the pulsar velocity dominates, an estimate for the 
pulsar transverse velocity based on the scintillation param- 



I^GHz Atd,s 



kms 



(13) 



where /vps is the ratio of the observer-screen dis- 

tance to the pulsar-screen distance and the p roportionality 
constant Ay = 3.85 x 10* (|Gupta et aLlligQ^ I. 

In Table [1] we give some basic parameters for the pulsar 
and observations as well as mean values for the basic scatter- 
ing parameters. The pulsar distance i s based on a parallax 
measurement by iBrisken et al.l l|2002f ) who also measured 
the pulsar proper motion, 19.5 ± 0.4 mas yr~^ which, with 
the distance, gives the transverse velocity listed in the table. 
The typical observation time for each dynamic spectrum and 
the total number of such observations are given. 

In the lower part of the table we give the mean values 
and rms scatter for the decorrelation bandwidth (Afd), scin- 
tillation time-scale (Aid) and the slope parameter (dt/du) 
derived from the ACF analysis. We note that the scintilla- 
tion timescale Atd is much less than the observation time (90 
min), so the effects of finite observation time discussed by 
(jStinebring et all |2000| ) should be minimal. The predicted 
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Figure 2. Contour plots of ACF for dynamic spectra shown in Fig. [T] There are 10 contours over the range zero to unity, successive 
contours are separated by an interval 0.1 and the dashed contours represent negative values. The observing MJD is shown in the top 
right corner of each panel. 



Table 1. Observation Parameters and Results 



Galactic longitude (deg) 


145.00 


Galactic latitude (deg) 


-1.22 


DM (pc cm-3) 


26.83 


Distance (kpc) 


1.06 ±0.12 


Vpm (km s-i) 


90 ±2 


Observation time (min) 


90 


Number of observations 


168 


(Az/d) (MHz) 


9.2 ±2.2 


(Aid) (min) 


17.1 ±3.3 


(dt/du) (min MHz-i) 


-0.36 ±0.46 


(h) 


8 ± 2 


C2 (10-4m-20/3) 


13.8 ±3.9 


u 


18.3 ±2.2 


Visa (km s~^) 


73 ± 17 



refractive timescale, AU based on the diffractive parame- 
ters (Equation [21), is about 4 d, but the observed value (Sec- 
tion l4.4|l is much less than that. The derived scattering level, 
Cn, and scattering strength, it, both indicate that scattering 
is strong along the path to PSR B0329+54. For the quoted 
scintillation velocity, a centrally located screen {x — 1) has 
been assumed. Our value is somewhat lower than the proper 



motion veloci ty. It is much lowe r than the value given for 
this pulsar bv lGupta et al.l l|l994h based on 408 MHz obser- 
vations, but they assumed a pulsar distance of 2.3 kpc. For 
the parallax distance of ~ 1.0 kpc, their velocity becomes 
105± 17 km s"'^, consistent with t he proper motion v elocity. 
From their 327 MHz observations iBhat etal] (llQQQal ) derive 
a scintillation velocity of 186 ± 17 km s~^ based on an as- 
sumed pulsar distance of 1.43 kpc. Scaled with the parallax 
distance, this becomes 155 ± 14 km s~^. While these results 
are not too discordant, it is clear that there are variations 
on timescales much longer than the nominal refractive time 
and also that the frequency scaling of the diffractive param- 
ete rs does not precis ely follow the Kolmogorov prediction 
(cf. IWang et al■ll2005^ . 



4 INTERPRETATION OF RESULTS 

4.1 Cross-Correlations Between Scintillation 
Parameters 

iRomani et al.l (|l986l ) suggested that correlations should ex- 
ist between variations of the flux density F, decorrelation 
bandwidth Az/j, and diffractive timescale Atd- These predic- 
tions assume a single-phase screen and a simple power-law 
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MJD-53076 

Figure 3. Time series of flux density (_F), decorrelation bandwidth (Ai^d), diffractive time-scale (Aid) and band drift rate (dt/dv) . The 
gap in the plots represents the interruption of the observations for about three days. 



Table 2. The theoretical predication and observational results of cross-correlation coefficients. 







{Ai/d, Aid} 


{Ai/d, F} 


{Aid, F} 


RNB86 


/3= 11/3 


0.75 


-0.76 


-0.50 




/3 = 4 


0.77 


-0.80 


-0.58 




13 = 4.3 


0.79 


-0.84 


-0.64 


SFM96 


T 


0.65 


-0.40 


-0.62 




Spearman Ts 


0.54 


-0.41 


-0.68 


BRG99b 


Spearman Ts 


0.04 


0.45 


-0.61 


This work 


Spearman Ts 


0.14 


0.30 


-0.07 




95% confldence interval 


-0.03 ~ 0.30 


0.17 ~ 0.44 


-0.22 ~ 0.08 



where Ri and Si are the ranks of the two quantities Xi and 
yi for which the correlation coefficient is computed and the 
summation is carried out over the total number of the data 
points. R and S represent the average values of Ri and Si. 
This method is less sensitive to the outlying points than 
the usual correlation coefficient r, which is computed with 
the same equation but with the actual data values rather 
than their rank, rs and its confidence interval are derived 
using the 'bootstrap' method as follows: assume we have an 
original object O with a set of N elements. We construct 
a new list with the same number of N elements from the 
original list by randomly picking elements from the list. Any 



description for the electron density fluctuations. In Table [2] 
we list the predicted correlati on coefficients for tu rbulence 
spe ctra of /3 =11/3, 4 and 4.3 l|Romani et~ai]|l986l) . 

IStinebring et al.l l|l996h and lBhat et"all l|l999bl ) studied 
the correlations between the variations of F, Ai^d and Aid 
for PSR B0329-I-54. Both papers use the Spearman rank- 
order correlation coefficient (r-s) to describe the correlation 
between the parameters, is given by: 



(14) 
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Figure 4. Scatter plots of Au^ vs At^, Av^ vs F and At^ vs F, respectively. 



one element from the list can be picked any number of times. 
The process is repeated 10000 times to give new series {Oi, 
O2....O10000} from which we obtain a statistical distribution 
of Vs. The mean value and 95% confidence limits are derived 
from this distribution. The res ults for Va are p r esente d in 
Table m along with results from lStinebring et al.l (|l996l ) and 

iBhat et al.l (Il999bh for comparis on. 

Although IStinebring et al.l (|l996l ) found significant 
cross-correlations between the scintillation parameters 
which were appro ximately in accord with the theoretical 
predictions, neither lBhat et al] (|l999lj ) nor our observations 
show such correlations. As Fig. |3]shows, little or no correla- 
tion is observed between Avd and Atd and the correlation 
between Avd and F even has the opposite sign to the the- 
oretical prediction. In contrast to the other studies, we ob- 
serve only a very weak correlation between Atd and F. The 
observed correlation s are clearly not we ll described by the 
thin-screen model of lRomani et al.l l|l986l) . Their predictions 
are not strongly dependent on the assumed slope of the fluc- 
tuation spectrum so other factors, for example, an extended 
scattering region, must be important. However, there have 
been no predictions for the expected correlations in this case. 



4.2 Structure function analysis 

We use a structure function analysis to study the refractive 
variations in the flux density and the diffractive timescale 
and decorrelation bandwidth. The structure function is de- 
fined to be: 



Table 3. Structure function noise estimates 



D{n) 



1 



■ ^ w{i)w{i + n) [F(i) - F\ 



i + n) 



(15) 

where D{n) is the structure function at lag of n units of the 
observation time (90 min), F{i) is the flux density value for 
the ith observation, 'u;(i) = 1 if F{i) exists and zero other- 
wise, {F) is the mean flux density, and N{n) is the number 
of products in D{n) for which w{i)w{i + n) 7^ 0. Typically, 
the structure function has three regimes: a noise regime at 
small lags, a structure regime characterized by a linear slope 
on a log-log plot and a saturation regime where it flattens 
out at large lags. The shape of the structure function may 
be used to derive the refractive modulation index and 
the refractive scintillation time-scale At^. The modulation 
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index rrir is expressed as mr = ^ D{oo)/2, where D{oo) is 
the saturation value of the structure function, and the re- 
fractive time-scale Atr is the time lag at which the structure 
function reaches half of its saturation value. The logarith- 
mic slope of the structure function, 7, is defined in the lin- 
ear structure region. Structure functions for the refractive 
variations in diffractive timescale and bandwidth may be 
similarly defined. 

Structure functions must be corrected for the effects 
of noise in the time series. The noise contribution to the 
structure function values, Dnoiso = ^a^^i^c = 2{a^st + crnicas)5 
where (Test is the statistical estimation error resulting from 
the finite number of scintles in the dynamic spectrum (Equa- 
tion |4]), CTmcas is the mean fractional uncertainty in the 
measured parameter. In Table [3] we list measured values of 

(Tmcas , (Test , <Tnoisc 

and Dnoisc for the three quantities, along 
with the uncorrected value of D{1), the 90-min lag value. 
Corrected values of the derived structure functions, that 
is, with Dnoise Subtracted, are shown in Fig. [5] Uncertain- 
ties on the structure fu nction values are given by a pin) = 
anoisc[8D(n)/N{n)Y^^ (|Stinebring fc CondonI [l990l ) . In all 
cases, the corrected value of -D(l) is below the extrapolation 
of the linear part of the structure function, indicating that 
we may have somewhat over-estimated the noise contribu- 
tion. 



4.3 Modulation Indices of Scintillation 
Parameters 

In strong scattering, the d iffractive para meters are modu- 
lated by refractive effects. iRomani et al. I (Ii986) gave rela- 
tions for the expected modulation indices for refractive vari- 
ations of flux density (rrir), diffractive bandwidth (rrib) and 
diffractive timescale (mt) for power- law fluctuation spec- 
tra with (3 = 11/3, 4 and 4.3. For the steeper spectra, 
the modulation indices are largely independent of scattering 
strength and frequency. They found that the depth of mod- 
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Table 4. Modulation indices for DISS parameters 
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Figure 5. Log- log plots of structure functions for observed vari- 
ations of: (A) Flux density _F, (B) decorrelation band-width Av^ 
and (C) diffractive time-scale Atj. The plotted structure func- 
tions have been corrected for the estimated noise contribution. 
The horizontal dashed lines are at the level corresponding to the 
observed noise-corrected modulation index. 



ulation is lowest for the Kolmogorov spectrum and larger 
for larger values of /3. Predicted modulation indices for PSR 
B0329- j-54 are given in columns 2 - 4 of Table [4] . Based 
on the lRomani et al.l ((1986,) theory, iBhat et all (jl^i) gave 
the modulation indices of F, Afd, and Afd for P = A and 
13 = 4.3. We list them in column 3 & 4 of Tabled 

The noise-corrected modulation indices may be com- 
puted directly from the time series using 



(16) 

where a is the observed rms fluctuation relative to the mean. 
The modulation index has an uncertainty given by 

1/2 
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0.31 ±0.03 
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0.16 ±0.02 



where A^scint is approximately the data span divided by the 
refractive timescale and A'obs is the number of independent 
observations. Using the noise estimates given in Table |31 the 
derived modulation indices and their uncertainties are given 
in column 5 of Table U The derived values are consistently 
lar ger than the predicted v alues for a Kolmogorov spectrum 
(cf. IStinebring et aLllioOol ). indicating a /3 ~ 4.0. 

Another possible intepretation of the relatively large ob- 
served modulation indices is that the scattering medium is 
extended so that the thin-screen approximation is not valid. 
The flux density modulation index for a statistically uniform 
scattering medium with a Kolmogorov spectrum {/3 = 11/3) 
co vering the whole path to the pulsar h a s been discussed 
bvlColes et al.l JigSThTCupta et all (|l993l '). IStinebring et all 
feOOd ') and lSmirnova et all (|l998i rwith the resuh 



1.05(Ai/d/i/) 



(18) 



For our observations of PSR B0329+54, this relation gives 
a predicted my ~ 0.45, larger than the observed value. This 
implies that, as an alternative to a steeper spectrum, the 
observed modulation index could be accounted for by an 
extended Kolmogorov scattering medium covering just part 
of the path. 

4.4 Refractive Timescales 

The horizontal dashed lines on the structure function plots 
(Fig. [5]) are at D{oo) = 2m^ , where m is the observed mod- 
ulation index given in in Table |4] The refractive timescale 
Ati is defined to be the lag at which the structure func- 
tion falls to 0.5D(oo). In the flux density case, this is at 
8 ± 2 h. The structure functions for the diffractive band- 
width and timescale are less well defined, but clearly also 
imply a short refractive timescale. The measured timescale 
is much less than the predicted value of about 96 h based 
on a central thin-screen model and the observed diffractive 
timescale (Table [l|. For an extended scattering medium 
with a Kolmog o rov s pectrum, evaluating the equations of 
ISmirnova et al.l (|l998l ) with practical units gives 



Atr.h = 278- 



R 



1/2 



A 



(19) 



(17) 



Substituting the parameters from our observations of PSR 
B0329+54 gives a value of about 32 hours, again much larger 
tha n the observed v a lue. 

ISmirnova et al.l (|l998t ) also consider the case of a thin 
screen at an arbitrary location along the path and show 
that Ati is proportional to the pulsar-screen distance (rps) 
whereas Aid is proportional to R/{R — Tps). Therefore, if 
rpa « R, refractive times are much reduced, whereas Aid is 
only weakly dependent on rps. This breaks the approximate 
proportionality of Atr and Atd (Equation [5} which applies 
to both a central thin screen a nd an extended medium. The 
scintillation arc observations of lPutnev fc Stinebringj (|2006l ) 
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suggest that most of the scattering for PSR B0329+54 oc- 
curs relatively close to the pulsar, supporting this interpre- 
tation of the short refractive timescale. 



4.5 Structure Function Slope 

The slope of the linear part of the structure function, be- 
tween the noise and saturation regimes, is related to the 
power-law index of the fluctuation spectru m. For a single 
thin screen, the slope should be 2.0 for /? < 4 H omani et al.l 
1 19861 : i mirnova et al.l Il998|). For an exte nded screen, the 
predicted slope is /3 — 3 ( Smirnova et al.l |1998 ) . Although 
the observed slopes are not very well determined, they are 
clearly <^ 2.0. The measured value for the flux density struc- 
ture function slope is 0.4 ± 0.2 and the values for the other 
two structure functions are similar. This is clearly inconsis- 
tent with scattering by a single thin screen. As we have ar- 
gued above, there are good arguments supporting the idea 
that the scattering screen in this direction is extended; in 
this case, the implied value of /3 = 3.4 ± 0.2, somewhat 
less than but marginally consistent with the Kolmogorov 
value. Based on the thin-screen model, observed diffractive 
bandwid ths and a r e fractive inter pretation of band slopes 
{dt/dv), iBhat et all l|l999al ) and Iwang et al.l (|2005l ) also 
obtain /3 values less than the Kolmogorov value for PSR 
B0329-I-54 and sever al other pulsars. Our re sult also agrees 
well with that from IStinebring et al.1 (|200G| ) who obtain a 
structure function slope of 0.5 ±0.1 from their flux density 
monitoring of this p ulsar. However, it contrasts with that of 
IShishov et al.l l|2003i ) where a structure function slope ~ 1.5 
is derived from observations at several different frequencies 
and interpreted in terms of weak plasma turbulence, again 
with a Kolmogorov spectrum. 



5 CONCLUSIONS 

Quasi-continuous observations of PSR B0329+54 at 1540 
MHz over 20 days have been analysed to investigate the re- 
fractive modulations of the pulsar flux density and diffrac- 
tive scintillation parameters. The data set was split into 
more than 150 individual observations, each of 90-min dura- 
tion, and two-dimensional auto-correlations and secondary 
spectra computed. Time series of the pulsar flux density 
(-F), decorrelation bandwidth (Afd), diffractive scintillation 
time-scale (Aid) and drift rate of features in the dynamic 
spectra (dt/dv) clearly show refractive modulation, although 
no evidence for arc structure is seen in the secondary spec- 
tra. 

Observed cross-correlations between variations in the 
flux density and diffractive parameters were much smaller 
than predictions based on the thin-screen model. In one case 
the sign of the correlation was opposite to the predicted 
value. Similar results have been obtained by other observers 
although different observations appear to give conflicting re- 
sults. 

In accordance with previous work, observed modulation 
indices are greater than predicted for a thin screen with a 
Kolmogorov fluctuation spectrum. This could be accounted 
for either by a steeper fluctuation spectrum (/3 ~ 4) or by 
scattering in an extended medium. The predicted modula- 
tion index for a scattering medium covering the whole path 



to the pulsar is greater than that observed, suggesting that 
the scattering medium, while extended, does not cover the 
whole path. Structure functions derived from the observa- 
tions indicate a short refractive timescale, 8± 2 h, much less 
than predicted from the thin screen model, and have a rela- 
tively flat slope, 0.4 ±0.2, again inconsistent with scattering 
by a thin screen. 

The observed modulation indices, structure function 
slopes and short refractive timescales all are consistent with 
scattering by an extended region with a Kolmogorov fluctua- 
tion spectrum which is concentrated toward the pulsar. This 
idea is supported by recent high-sensitiv i ty ob servations of 
PSR B0329±54 bv lPutnev fc Stinebrin3 l|2006l ) which show 
indistinct scintillation arcs corresponding to extended scat- 
tering regions relatively close to the pulsar. 
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